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ABSTRACT 

We report on the merger-induced generation of a shock-heated gas wind and formation of a remnant gas halo 
in simulations of colliding disk galaxies. The simulations use cosmologically motivated initial conditions and 
include the effects of radiative cooling, star formation, stellar feedback and the non-adiabatic heating of gas. 
The non-adiabatic heating, i.e. shocks, generated in the final merger forces gas out of the central region of the 
merger remnant and into the dark-matter halo. We demonstrate that the amount of heating depends on the size of 
the progenitor disk galaxy as well as the initial orbit the galaxies are placed on. Based upon these dependencies, 
we motivate a possible recipe for including this effect in semi-analytic models of galaxy formation. 
Subject headings: galaxies:evolution — formation — interactions — ISM — methods:N-body simulations 



1. INTRODUCTION 

Galactic winds are a common phenome non associated 
wi th starbursting galaxies in the loca l ( Martin Il999t 

f ickman et alJ2000tlStrickland et alJ2002l) and high redshift 
;ttini et al. 2002; Adelberger et al. 2003) universe. It is gen- 
erally assumed that feedback from massive stars, either stellar 
winds or supernova ejecta or both, power these outflows of hot 
gas. The starburst origin of galactic winds in nearby galax- 
ies is supported by the observed correlation bet ween wind- 
induc ed mass loss and the star fo rmation rate ( M artiri 119991: 
Rupk eTVeilleux. & Sandersl20 02). In addition to transporting 
mass and energy from the centra l star-forming region, galac- 
tic winds also transport metals (Heckman et al. 2000). De- 
pending on the energy in the winds and the depth of the halo 
potential well, outflows such as these may be responsible for 
heating and polluting the intergalactic medium. 

Galactic winds likely play a significant role in galaxy for- 
mation and evolution. As a source of energy, or feedback, 
galactic winds regulate the conversion of gas into stars. In 
recent semi-analytic models of galaxy formation, winds (in 
addition to a UV background) are required to remove gas and 
suppress star formation in small mass halos, making the theo- 
retical halo mass func tion consistent with the observed galaxy 
luminosity function ISomervillel 120021 iBenson et all 2003). 
Energy input due to galactic winds counteracts the efficient 
gas cooling in galactic halos - the so-called "overcooling" 
problem of theoretical models or simulations. Interestingly, 
not only do these models produce too many faint galaxies, but 
they typically produce too many bright galaxies as well. Thus, 
overcooling is a problem in both low and high mass halos. 

While feedback from massive stars provides a physically 
reasonable and effective mechanism by which the low mass 
end of the halo mass function can be reconciled with the 
galaxy luminosity function, it is unclear that supernova-driven 
winds can provide enoug h energy to s o lve th e overcooling 
problem for large halos. Benson et al. (2003), using semi- 
analytic models, were able to match the bright end of the lu- 
minosity function only by assuming extremely efficient ther- 
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mal conduction or by including galactic winds whose energy 
greatly exceeded that produced by supernovae. 

It is clear additional sources of energy are required during 
the process of galaxy formation and in this letter we suggest a 
possible candidate. We demonstrate that simulations of disk- 
galaxy mergers generate shock-heated winds which can put 
gas with thermal energy up to 1 keV per baryon into the dark- 
matter halo. These winds are not induced by supernovae, but 
instead caused by orbital energy which is transferred to the 
diffuse gas. Since major mergers are more common during the 
hierarchical formation of massive structures, energy input due 
to this process will preferentially act in more massive halos. 
Further, the generation of hot gas due to the merging process 
suggests that we should be able to detect hot gas in nearby 
merging systems. 

The structure of this letter is as follows: ^summarizes our 
simulation techniques and our initial conditions, [j3l presents 
an overview of the galaxy merger process focusing on the con- 
ditions relevant to generating hot galaxy winds, Ogives prop- 
erties of the outflow material and presents a physically moti- 
vated fitting formula to capture our results which we discuss 
in 

2. SIMULATIONS 

We use GADGET iSpringel. Yoshida. & White! 1200(1 
to simulate the merger between two identical copies 
of our model disk galax y. Based on the work of 
Springel & Hernauisn ( 120021) . and adopting their terminology, 
we use the the conservative-entropy version of smoothed par- 
ticle hydrodynamics (SPH), although we remark that because 
we include shocks and radiative cooling entropy is not con- 
stant throughout the simulation. Gas is assumed to be a pri- 
mordial plasma and can cool, form stars and produce super- 
novae which release energy and return metals to the ISM. Star 
formation is based upon the local dyna mical time scale. Stel- 
lar feedback is implemented similarly to Springel (2000); we 
define a feedback energy reservoir in which supernova feed- 
back energy is stored. This reservoir has a dissipative time 
scale longer than the cooling time and can provide pressure 
support in star-forming regions. Temperature and pressure are 
calculated as the sum of the thermal and feedback reservoirs. 
We note that we have varied our star formation and feedback 
assumptions and the results we report here ar e robust. Details 
of our methods are given in lCox et"ail (|2004). 

Star-formation and feedback parameters are adjusted so that 
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Table 1. Disk Galaxy Initial Conditions 



Label 


Mvir 


Rvir 


c 


A 




/ 


R d 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


5xMW 


700 


499 


11 


0.05 


0.04 


0.5 


10.6 


MW 


140 


291 


12 


0.031 


0.04 


0.1 


3.4 


Dwarfl 


14 


135 


15 


0.05 


0.04 


0.5 


2.1 


Dwarf2 


1.4 


63 


20 


0.05 


0.05 


0.5 


1.0 



NOTE. — Column 1 : Disk galaxy label. Column 2: Galaxy virial mass in units of 
10 10 Mq ; Column 3: Virial Radius R v ir, in kpc; Column 4: Initial halo concentration 
c=R v j r /r s where r s is the NFW scale radius; Column 5: Dimensionless spin parameter; 
Column 6: Fraction of virial mass in an exponential baryonic disk; Column 7: Fraction 
of exponential baryonic disk mass in a collisional gaseous component; Column 8: Disk 
scale length R^, in kpc. 



an isolated disk gala xy's star f ormation follows the empirical 
Schmidt law (Kennicutt 1998). Once these parameters have 
been set, we collide two disk galaxies by fixing them on a 
prograde orbit defined by the eccentricity e and the pericen- 
tric distance i? pcr i, as determined from the initial energy and 
angular momentum. One disk is in the orbital plane while the 
other is inclined by 30°. 

We simulate mergers between four different identical pairs 
of progenitor disk galaxies. The largest represents a disk of 
roughly 5x the mass of the Milky Way, which is our sec- 
ond largest progenitor galaxy. We select parameters for the 
Milky Way from model Al of Klvpin .~Zhao. & Somervillel 
(2002) and for simplicity ignore the bulge component. Our 
final two progenitors are dwarf galaxies of one-tenth and one- 
hundredth the mass of the Milky Way. Initial disk galaxy pa- 
rameters are listed in Tabled 

Each galaxy model consists of an exponential disk, of scale 
length Rd and gas fraction /, wh ich resides in the center 
of a spherically symmetric NFW (Na varro. Frenk. & White! 
1997) dark matter hal o. The nu merical construction of this 
system is described in Springel (2000). We use 180,000 to- 
tal (40,000 gas) particles for our largest progenitor galaxy, 
70,000 (20,000) for our Milky Way-like galaxy and 30,000 
(5,000) for each of the dwarfs. GADGET smoothing lengths 
are 570 pc for dark matter particles and 140 pc for bary- 
onic particles in our two larger systems while the dwarfs have 
smoothing lengths of 140 pc and 57 pc. Identical copies of 
each galaxy model (as specified in Tabled are put on a vari- 
ety of orbits which are listed in Table|2] Each interaction was 
simulated until the merger was relaxed, typically 700 Myr af- 
ter the nuclei coalesce, and total simulation times were be- 
tween 2-5 Gyr. 

3. GALAXY MERGERS 

The morphology and induced star formation of gaseous 
g al ax y collisions has been discussed in a number of pa- 
pers (iBarnes & Hernauistl Il996t iMihos & Hernauistl Il996t 
Springe] 2000). Here we briefly emphasize the aspects of the 
merger process relevant for our purpose. The simulation be- 
gins with a wide separation between the two galaxies and they 
are essentially unaffected by each other's presence. As the 
galaxies reach pericenter the disks become tidally distorted 
and a small burst of star formation ensues. At this point dy- 
namical friction has extracted orbital energy from the origi- 
nally unbound orbit but the galaxies are still able to briefly 
separate before they finally merge. At the final merger the 
two disks become effectively destroyed and we are left with a 
spheroidal-looking object. 
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FIG. 1 . — Velocity field of the gas overlaid on a gray scale map of the gas 
density in a 7 kpc slice through a plane perpendicular to the orbital plane at 
2.4 Gyr after the start of the simulation and 200 Myr before the final merger. 
The length of the arrow in the upper left corner represents 200 km/s. From 
this view the coherent wind motion is clearly visible as it escapes the high- 
density merging region out of the orbital plane, i.e. moving along the path of 
least resistance. 



We find that the total fraction of gas converted to stars dur- 
ing the simulations rang ed from 68% to 92% in line with 
previous studies (Mihos & Hernquist 1996; Springel 2000). 
While the star formation is generally efficient at converting 
gas to stars, we note that the bursts of star formation dur- 
in g the final merge r are less efficient than in previous work. 
In C ox et alJ §004) we show that this results from using the 
conservative-entropy version of SPH. Each orbit has a unique 
star formation history which dictates how much star formation 
occurs quiescently or during tidal triggered episodes, but we 
note that the total gas consumption is only weakly correlated 
with orbit. 

At the time of the final merger between the disk galaxies, 
the majority of remaining gas is cold and of relatively low 
density. Most of this gas was in the outer portions of the orig- 
inal disks. Some of this gas was ejected in long tidal tails 
during the first encounter of the two disk galaxies, but by the 
final merger much has fallen back to its parent galaxy. As the 
two primarily cold gas disks attempt to follow their dominant 
dark halos into the final merger, the collisional nature of the 
gas begins to dominate and a large amount of gas shocking oc- 
curs. These shocks can heat the gas to temperatures between 
10 6 and 10 7 K. 

The hot gas then expands adiabatically and finds little re- 
sistance out of the orbital plane as shown in Figure ^ Gas 
collectively moves at speeds > 200 km/sec, and quickly fills 
the dark-matter halo. It is often the case that there are several 
episodes of gas ejection similar to Figure^due to the multiple 
shock events as the nuclei oscillate and finally merge. Since 
we have not included a gaseous halo in our initial conditions 
we do not capture the interaction of the wind with the ambient 
gas within the original halo. 

Within the gaseous component, the final merger corre- 
sponds to a transfer of bulk kinetic energy to thermal en- 
ergy, indicative of shocks. While the amount of energy gener- 
ated by supernovae associated with star formation during the 
course of the merger is quite large, its contribution is local- 
ized to the high density regions and spread out over time. The 
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Table 2. Galaxy Collision Orbits 
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Label R sta rt (Rperi, e) 



5xMW 429 (1.7, 0.99), (7.1, 1), (21.4, 1), (35.7, 0.9), (35.7, 1), (85.7, 0.6) 

MW 286 (0.1,0.99), (7.1, 1), (21.4, 0.9), (35.7, 0.8), (178.6,0) 

Dwarfl 143 (1.1, 1), (2.9, 1), (5.7, 1), (10, 1), (15.7, 1), (21.4, 1), (31.4, 1) 

Dwarf2 57 (0.6, 1), (1.1, 1), (2.4, 1), (4.1, 1), (6.4, 1), (8.9, 1), (12.9, 1) 



NOTE. — Column 1: Disk galaxy label. Column 2: Initial galaxy separation R s tart> ln kpc; Column 3: The pericentric distance R pcr i, in kpc, and the eccentricity of each orbit 
simulated. 



supernovae feedback energy contributes a negligible amount 
to the energy and entropy gain in the lower density gas at the 
final merger. During the final merger the majority of the en- 
tropy gain is due to low-density gas, which is shock heated 
while never (or only briefly) reaching star-forming densities. 
This emphasizes that the winds reported here are primarily 
due to shocks and not star formation. One might expect the 
thermal energy to be efficiently radiated away, but in practice 
the gas quickly accelerates out of the central region to densi- 
ties where the cooling time ranges from > 300 Myr to several 
Gyr. 

We verified that resolution effects were not limiting our 
ability to capture the shock heating during the galaxy mergers 
by resimulating the MW(R per i=7.4,e=l) orbit with a factor of 
10 higher and lower mass resolution. The amount of hot gas 
in the merger remnant and the average energy of this gas var- 
ied by less than 3% in these resimulations, emphasizing that 
we are not limited by resolution. We speculate that cosmolog- 
ical simulations do not see shock heating in mergers because 
they efficiently form stars from all cold gas inside dark matter 
halos long before the dark matter halos finally merge. This 
will be investigated in a fut ure paper. 

The "phase difference" dNavarro & Whitelll993l) between 
the dominant collisionless center of mass and the gas during 
the final merger is dictated by the orbits the initial galaxies are 
placed on. Nearly radial orbits tend to have their nuclei oscil- 
late through each other during the final merger providing for 
ample phase lags and hence ample shock heating. Conversely, 
circular orbits demonstrate very little gas heating during the fi- 
nal merger because the nuclei spiral around each other as the 
galaxy centers coalesce. In this case gas from each disk quies- 
cently 'slides' around the other disk's gas and little shocking 
takes place. 

4. HOT GAS IN MERGER REMNANTS 

In Figure [2] we detail the thermal energy per proton of 
all gas in each merger remnant. There is a trend for orbits 
with smaller impact parameters, here called b and defined as 
Rpcri/Rd, to generate gas with more thermal energy. Since 
our initial conditions did not provide any low density medium 
representing ambient gas we are not capturing the inevitable 
process by which the out-flowing gas can thermalize any of 
its kinetic energy. We do not show this kinetic energy in Fig- 
ure |3 since there are additional uncertainties with regard to 
energy loss due to radiation or the production of cosmic rays, 
but we remark that the wind's kinetic energy is approximately 
10-50% of its thermal energy. 

From energy considerations alone we can attain an approx- 
imation of the gas thermal energy plotted in Figure [2] While 
the initial kinetic energy and total orbital energy are depen- 
dent on the assumed eccentricity and pericentric distance, the 




FIG. 2. — Remnant gas energy in keV per proton as measured 500 Myr 
after the merger between two disk galaxies plotted as a function of the or- 
bital pericentric distance R per i divided by the original disk scale length R^. 
Open circles are our the progenitor disk galaxy 5xMW, the largest progenitor 
disk galaxy, solid squares are MW disks, the asterisk are Dwarfl and solid 
triangles are Dwarf2. Solid lines are fits to the data using equation and 
described in the text. 



initial potential energy per unit mass is 



GM vir 

Rstart, 



0.22 keV /proton 



\Q 12 Mr. 



/200fcpc\ 

V Rstart / 



independent of the orbit. During the merger this energy is 
absorbed by the galaxies, consistent with Figure[2] 

For a more complete understanding of the origin of Figure|2] 
we must consider how shocks thermalize kinetic energy. The 
initial disk gas is cold (temperature ~ 10 4 K) prior to shock 
heating and the sound speed is < 15 km/sec. The orbits we 
assume dictate the typical velocity with which the galaxies in- 
teract to be > 100 km/sec and hence the gas in each progenitor 
disk will undergo strong shocks (Mach number > 10) as the 
disks collide during the final merger. The Rankine-Hugoniot 
jump conditions in the strong shock regime give 



itial' 



(2) 



where Tfi na i is the temperature of the gas after the shock, 
and mvf nitial is the initial kinetic energy of gas. In physi- 
cal terms the jump condition states that the initial kinetic en- 
ergy is transfered to gas thermal energy via the strong shock 
with almost 50% efficiency. Moreover, the strong shock only 
occurs when the progenitor disk galaxies pass through each 
other rather than orbit around each other, and thus we reason 
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it is the initial kinetic energy in the radial direction which gets 
converted to heat because of shocks, i.e. Vimtiai ~ Uradiai- 

It is straightforward to determine the initial kinetic energy 
along the line connecting the centers of the initial disk galax- 
ies as 



v 2 

^radial 



GM V 



1 



(1 



R 



(3) 



pen 



where R pcr i and e are the pericentric distance and orbit eccen- 
tricity, respectively. We see that the initial radial kinetic en- 
ergy is inversely proportional to the pericentric distance, and 
hence initial orbits with large impact parameters should have 
much less gas shocking and thus less hot gas in their merger 
remnant, as seen in Figure[2] 

Based upon the above assumptions we propose a model by 
which we can capture the general features of generating hot 
gas by disk galaxy major mergers. Primarily the generation 
of hot gas is inversely proportional to the impact parameter 
b defined above. We find that the following simple formula 
accurately captures this relationship, 



fcTfi na i/proton 



^mcrg 



b + R d /R A • 



(4) 



where e mer g is a constant with units energy per proton, b is the 
impact parameter defined above, Rd is the disk radial scale 
length and R^ is a distance which acts as a fudge factor to 
scale the heating dependence on the disk size. This formula 
is fit to the data using a least squares fitting routine and dis- 
played in Figure |3 For mergers of the 5xMW and MW pro- 
genitor disk galaxies we find e merg scales with the mass as 
expected from equation (0, its value can be approximated as 



IkeV/proton 



(5) 



and Ra is roughly 1.18 kpc. Fits for the smallest progeni- 
tors, Dwarf 1 and Dwarf2, are not well approximated by this 
formula as the thermal energy of all remnants was approxi- 
mately 3 eV/proton and 1 eV/proton, respectively. We expect 
a fair amount of scatter due to varying gas distributions and 
gas mass fractions in the progenitor disk galaxies. More sim- 
ulations must be performed over a wide range of initial con- 
ditions to determine if equation (0) or Q better captures the 
full thermalization of orbital energy due to the merger. 

Very little hot gas is energetic enough to escape the poten- 
tial well of the halo, and the majority comes into equilibrium 
with a profile reasonably well repres ented by a spherical (3- 
profile Eke. Navarro. & Frenkl ( 119981) . The entropy per pro- 



ton generated via the merger process is < 400 keV cm 2 and 
follows a trend very similar to that seen in Figure|2] 

5. DISCUSSION 

Hydrodynamic simulations of colliding disk galaxies 
demonstrate the ability to generate a significant amount of 
hot gas from shocks which abundantly occur in the central 
several kpc region during galaxy collisions. As the centers of 
disk galaxies coalesce, collisional gas cannot inter-penetrate 
and gets heated while attempting to follow the potential well 
dominated by the collisionless dark matter and stellar com- 
ponents. We note that any isother mal gas assumption used 
to simplify the treatment of gas (Barnes & Hernauist 1996; 
Mihos & Hernquisl f 1 996 ) would not reproduce the process 
we discuss here. 

In general, hot gas generated by galaxy major mergers 
could be an important source of energy during the process 
of galaxy formation. Not only do radial orbits pump energy 
and entropy into the gas but they provide an effective mech- 
anism to recycle moderately dense cold gas into diffuse hot 
gas, providing a feedback loop for energy, metals and mass 
that correlates with the hierarchical build-up of any galaxy. 
In reality the process we describe here may be augmented by 
contributions from stellar winds, supernovae, and AGN which 
will increase the total gas heating. The remnants we report 
here have cooling times of several Gyr and could provide an 
enriched medium with which a gas disk could subsequently 
be formed. 

Our results suggest that any galaxy system which has been 
involved in a major merger should have hot gas residing 
throughout its dark matter halo. Additionally, the tempera- 
ture of this gas, at a fixed time after the final merger, should 
correlate with the mass of the progenitor disk galaxies and 
their merger orbit. Further work will investigate if the pro- 
files in our merger remnants resemble the extended gas dis- 
tributions fojmd_Jn_jekxed_ejhjticjddominated groups as 
found by Mulchaev & Zabludoff ( 1998) or if the outflowing 
gas is similar to the diffuse gas found in ongoing mergers 
such as Arp220 (McDowell et al. 2003) and 'The Antennae' 
( Fabbian oeTalteOOll) . 
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